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About 40-50 % of the quarkonium ground states J/ip(lS) and T(1S) produced in 
hadronic collisions originate from the decay of higher excitations. In a hot medium, these 
higher states are dissociated at lower temperatures than the more tightly bound ground 
states, leading to a sequential suppression pattern. Using new finite temperature lattice 
results, we specify the in-medium potential between heavy quarks and determine the 
dissociation points of different quarkonium states. On the basis of recent CDF data on 
bottomonium production, we then obtain first predictions for sequential T suppression in 
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nuclear collisions. 
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1. Introduction 



The large values of the charm and bottom quark masses permit potential theory to 
provide a realistic account of quarkonium spectroscopy []|]-0- From these studies, it is 
known that the intrinsic length scales of quarkonia are much smaller than those of the 
usual hadrons, with rj/^ ~ 0.2 fm and ry ~ 0.1 fm for the radii of the lowest cc and 
bb vector mesons, respectively, in contrast to about 1 fm for light qq state radii. This 
"short-distance" nature of quarkonium states suggests that at least some features of their 
production in hadronic collisions should be accessible to perturbative QCD calculations, 
and that indeed turns out to be the case. 

The simplest and most general model for quarkonium production, the color evapora- 
tion model [Q, postulates that the cross-section for the production of a given charmonium 
or bottonium state is simply a fixed (energy-independent) fraction of the corresponding 
perturbatively calculated cc or bb production cross-section. The resulting predictions for 
the energy variation of the J/ip and T hadroproduction cross-sections are experimentally 
very well confirmed, for the T over a range from yfs ~ 20 to 1800 GeV ||. The assumed 
energy independence of the production ratios of 2S/1S and 3S/1S states is found to hold 
over the same range |J. 

Higher quarkonium excitations decay into lower states with generally known branching 
ratios and widths. As an example, the ip' = ip(2S) decays into J/ip(lS) + anything with 
a branching ratio of 55 %, after a mean life-time of more than 10 3 fm. As a consequence, 
J/i/j or T production in hadronic collisions occurs in part through the production of higher 



1 



excited states which subsequently decay into the quarkonium ground states. It is known 
experimentally that for both J/ip and T about 40 - 50 % of the hadroproduction rate is 
due to such feed-down from higher excitations [fj - []l(J . 

Quarkonium production through feed-down becomes particularly interesting when 
quarkonium states are used to probe the hot and dense medium created in high energy 
nuclear collisions. It was predicted that color deconfmement (quark-gluon plasma for- 
mation) would lead to J/ip suppression, since sufficiently hot deconfined media dissolve 



any cc binding JLT]. However, different quarkonium excitations will dissolve at differ- 
ent temperatures of the medium [|12| , and through a lowering of the open charm (beauty) 
threshold with temperature, dissociation by DD or BB decay becomes possible for higher 
excited states even below the deconfmement point |I3| . Since the life-time of the exci- 
tations is much larger than that of the medium, feed-down production will result in a 
characteristic sequential suppression pattern [|TJ|, |T5| , with the fraction of J/ip or T pro- 
duced through the decay of higher excitations becoming suppressed at lower temperatures 
than the directly produced ground states. 

To fully predict this sequential suppression, two prerequisites are needed. We have 
to know what fractions of the quarkonium ground state production originate from which 
higher excitations, and we have to know at what temperature or energy density of the hot 
medium a given excitation dissolves. The first can be determined either experimentally 
or through a viable model for quarkonium hadroproduction. The second is a well-defined 
problem for finite temperature lattice QCD studies. At present, neither problem is com- 
pletely solved. However, in the case of charmonium, the production rates of the higher 
excitations are experimentally known, and for bottonium production, recent Fermilab 
data |TO, Ofl provide the basis for fairly reliable estimates. This determines the structure 



(i.e., the sequence and the different heights) of the various suppression steps for J/ip and 
T production, but not the actual positions of these steps as function of the temperature 
or energy density. To estimate these, we make use of recent lattice studies calculating the 
temperature behavior of the heavy quark potential in full QCD |l(J. As a result, we 
obtain a modified J/ip suppression pattern [[14], |I5| , including the decay of ip' and Xc in 



confined matter |L3j , and then first quantitative predictions for the sequential T suppres- 
sion to be studied in forthcoming RHIC and LHC experiments. 

The structure of this paper is the following. In section 2, we first summarize the 
hadroproduction cross sections for the different charmonium states and then determine 
the corresponding bottonium cross sections from the mentioned new CDF data. As a 
result, we can fully describe the origin of the J/ip and T produced in hadronic collisions 
in terms of feed-down from higher excited states. In section 3, we consider the temperature 
dependence of the heavy quark potential obtained in recent lattice QCD studies. Solving 
the corresponding Schrodinger equation, we then determine in section 4 the dissociation 
parameters for the different states in a deconfined medium. This leads to quantitative 
estimates of the sequential J/ip and T suppression patterns. 

2. Quarkonium Production and Feed-Down 

A. Charmonium States 

It is well known that J / ip production in hadron-hadron collisions is to a considerable 
extent due to the production and subsequent decay of higher excited cc states f|-||. We 
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shall here summarize the situation following systematic studies using pion and proton 
beams at 300 GeV incident energy ||. In the first two columns of Table 1, we list the 
cross-sections af obtained for the direct production (excluding feed-down) of the different 
charmonium states ip(lS), x(l-P) an d ip'=ip(2S) in 7r~-nucleon and p-nucleon interactions, 
normalized to the overall measured J/ip cross-section aj/^, which includes all feed-down 
contributions. Hence Ri(ir~N) = (rf(7r~N)/aj/^(7r~N) for the directly produced state i 
in 7i~N interactions, and similarly for p N collisions. 

Making use of the branching ratios B[xx{lP) -> ip(lS)] = 0.27 ± 0.02, B\x 2 (lP) -> 
il>(lS)] = 0.14±.01, and B[ip(2S) -> ^(1^)] = 0.55±0.05, one obtains the fractional feed- 
down contributions of the different charmonium states to the observed J/ijj production; 
these are shown in the next two columns of Table 1. Also listed are the dissociation 
energies E l dis , 

E Ais = 2M D -M i . (1) 

measuring how far the mass of state % lies below the zero-temperature open charm thresh- 
old 2M D =3.740 GeV. 



state 




Riip N) 


fi{*-N) [%] 


fi(p N) [%} 


E dis [MeV] 


J MIS) 


0.57 ± 0.03 


0.62 ± 0.04 


57 ± 3 


62 ± 4 


0.642 




0.72 ± 0.18 


0.60 ± 0.15 


20 ± 5 


16 ± 4 


0.229 


X2(IP) 


1.04 ± 0.29 


0.99 ± 0.29 


15 ± 4 


14 ± 4 


0.183 




0.14 ± 0.04 


0.14 ± 0.04 


8 ± 2 


8 ± 2 


0.054 




1 


1 


100 


100 





Table 1: Cross-sections for direct charmonium production in rc~N and pN collisions, 
normalized to the overall J/ip production cross section in the corresponding reaction 0; 
feed-down fractions and mass gap to the open charm threshold. 

From Table 1 it is seen that some 60 % of the observed J/ip are directly produced, 
about 30 % come from x and about 10 % from ip' decay. According to the color evaporation 
model, feed-down fractions as well as cross section ratios are energy-independent. The 
results shown in Table 1 for the ratio of ip' to overall J/ip production are in excellent 
agreement with a variety of experimental results over the range = 18 — 65 GeV, which 
give 0.14 ± 0.034 as average ||. Even Tevatron results at y^i = 1.8 TeV for charmonium 



3 



transverse momenta pt > 5 GeV lead to 0.19 ± 0.05 and are thus within errors in accord 
with the quoted average value. Similarly, the ratio of x to J/i/; production is found 
to be constant over the range y/s = 7 — 65 GeV f5|. As noted, the color evaporation 
model postulates that the different charmonium state cross-sections af are constant (i.e., 
energy- independent) fractions c$ of the overall 'hidden charm' cross-section 

af(s) = aa c - c (s) with M C - C <2M D . (2) 

This assumption is thus indeed very well satisfied, even though the corresponding pro- 
duction cross-sections themselves vary considerably in the energy range in question, in 
accord with the perturbatively calculable variation of <r C c(s). Note that J2i Q << 1> since 
more than half of the cc pairs formed with M c5 < 2Me> acquire the missing energy from 
the color field and then contribute to open charm production. 

The color evaporation model does not predict the values of the factors q. More de- 
tailed arguments do, however, lead to some relations between the different cross-sections. 
Projecting a color singlet cc state onto different quantum number configurations leads to 



the estimate 18 



a d {lS) ~ T{tfj{lS) -> e+e-) \ My , 



0.24, (3) 



where T denotes the corresponding dilepton decay width of the state in question, M its 
mass. The values 0.23 ± 0.06 and 0.22 ± 0.05 obtained from Table 1 for n~N and pN 
collisions, respectively, agree well with relation (|3|). 

The ratios between the different XiO-P) states are predicted to be governed essentially 



by the orbital angular momentum degeneracy [[Lj|; we thus expect for the corresponding 
cross-sections 

Xo(lP) : Xi(lP) : XsO-P) = 1:3:5. (4) 

From Table 1 we have for ir~N collisions X2O-P) / XiO-P) — l-44±0.38 and thus reasonable 
agreement with the predicted ratio 1.67. For pN interactions, the experiment measures 
only the combined effect of Xi an d X2 decay (30 % of the overall J/ip production); the 
listed values are obtained by distributing this in the ratio 3:5. 

B. Bottonium States 

Production cross-sections for the different (nS) bottonium below the open beauty 
threshold are known over a considerable energy range ||, 0], and as for charmonia, the 
resulting ratios are with very good precision energy-independent || . Corresponding cross- 
sections for the sub-threshold (nP) states are so far available only for transverse momenta 
p T > 8 GeV/c fllP] . To analyse the complete feed-down pattern, we thus have to find a 
way to extrapolate these data to pt = 0. 

In Fig. |I|, we show the pr distributions of the T, T' and T" states. It is evident that 
they exhibit a very similar transverse momentum behavior, which can be parametrized 
quite well by the form 

da 1 

^- = iV J exp{(-0.415K} ) (5) 
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where iVj specifies the overall cross-section values of the different states i. We shall 
therefore now assume that all bottonium states, nP as well as nS, are governed by the 
same distribution, so if we know a production ratio in one pt inteval, we know the 
ratio of the overall cross-sections. To test the validity of assuming such a universal pr 
dependence, we have compared data for the different nS states taken in the interval 
1 < Pt < 18 GeV to an extrapolation of corresponding data for 8 < px < 18 GeV. The 
resulting difference of about 10 % indicates the uncertainty inherent in our procedure. 

To complete the feed-down analysis, we need the decay branching ratios from higher 
to lower bottonium states, which can be readily computed from the Particle Data Com- 
pilation. From the measured cross-sections including feed-down effects (denoted by a) 
we now want to reconstruct the direct production cross-sections (denoted by af) of all 
sub-threshold bottonium states i. This leads to the relations 

a(T)=J2B[i^T]af, (6) 

i 

o{T)=Y J B[i^V]at (7) 

i 

and 

a(T") = a d (3S) (8) 
for the three different bb S states. The branching ratios B[i — > T] and B[i — > T'] are 



compiled in [21 



The experiment [IQJ further provides the fractions F\ and F-f of T production 



coming from the different IP and 2P Xb states. For these fractions we know 

F^ p a(T) = £ B[i - (1P)]B[(1P) - T] af (9) 

i 

and 

Fy P a(T) = £ 5[i -> (2P)]B[(2P) -> T] af (10) 

i 

in terms of the direct production cross sections for the different states. The branching 



ratios are again given in pT| . Eqs. (H) and (|10| ) illustrate that the observed \b states 
which can decay into T arise themselves in part through feed-down from still higher 
excited states. We have to specify the different fractions for this, since e.g. the melting of 
the T(2S) will also remove that fraction of the observed Xb(lS) production which comes 
from T(2S) decay. 

The experimental values for F^ p and F^ p are given in [0 for transverse momenta 
Pt > 8 GeV/c. Since the different S states show a universal pt dependence, we assume 
the same to hold for the P states and thus take the measured values 

F\ p = 0.27 ±0.11 , Ff = 0.11 ±0.06 (11) 

to remain applicable for the entire pt range. Making use of the measured overall cross 
sections at = 26.9 ± 0.6 mb, ay = 13.1 ± 0.5 mb and ax" = 5.5 ± 0.5 mb, we can then 
solve the five equations ([]) - (10) for the different direct production cross sections. The 
results are shown in Table 2, normalized to a-f, together with the relative fractions each 
state contributes to the overall T production. For the P-states Xb, the shown values are 
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based on the overall cross-section a(xb) = c(x&o) +&{Xbi) + ° r (Xb2), with the three orbital 
states assumed to contribute in the ratios 1:3:5 |l9fl . Again we also list the values 

E dis = 2M B - Mi (12) 

of the corresponding zero-temperature dissociation energies. 



state 


Ri{pp) 


fi(pp) [%} 


E dis [GeV] 


fi (PP) NRQCD [%] 


T(1S) 


0.52 ± 0.09 


52 ± 9 


1.098 


0.52 ± 34 


Xb(lP) 

/ v, v \ J 


1.08 ± 0.36 


26 ± 7 


0.670 


0.24 ±8 


T(2S) 


0.33 ± 0.10 


10 ± 3 


0.535 


8±7 




0.84 ± 0.4 


10 ± 7 


0.305 


14 ±4 


T(3S) 


0.20 ± 0.04 


2 ± 0.5 


0.203 


2±2 


T 


1 


100 




100 



Table 2: Cross-sections for direct bottomonium production in p — p collisions, normalized 



to the overall T production cross section (TO, B[; feed-down fractions and mass gap to the 



open bottom threshold; feed-down fractions obtained in NRQCD. 

For the feed-down pattern in T production we thus find approximately 50 % direct 
T(IS') production, 30 % from direct Xb(l-P), 10 % from direct T'(2S) and 10 % from direct 
Xb(2P) production and subsequent decay into T(IS'). Again the energy independence 
of these fractions, as required by the color evaporation model, is well satisfied for the 
measured ratios 0. To test the consistency of the feed-down pattern with more recent 
theoretical considerations, we have also calculated the fractions for the whole measured 



Pt interval using the NRQCD factorization formula |2C, 21 1. The results are presented in 



the last column of Table 2; the details of the calculations are given in the appendix. 

3. The Heavy Quark Potential in Hot Media 

In finite temperature lattice QCD, the temperature behavior of the static QQ potential 
V(T,r) is obtained from Polyakov loop correlations measuring the free energy F(T,r), 

- T In < L(0)L+(r) >= F(T, r) + C = V(T, r)-TS + C, (13) 
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where S denotes the entropy due to the introduction of an unbound cc or bb pair into the 
medium and C the (undetermined) Polyakov loop normalization. This can in principle 
be fixed by requiring that at very short distances, r « T -1 , the potential has the purely 
Coulombic form a/r, since in the limit r — > 0, the effects of the medium should become 
negligible. At present, however, lattice calculations for high temperatures are probably 
not yet precise enough to reach the small r range required; for T < T c , the normalisation 



is found to be more reliable |13|]. We therefore here leave C open; as will be seen, this 
does not affect the determination of the dissociation points of the different bound states; 
it only prevents a reliable determination of the binding energy. The definition of S in Eq. 
flOp assumes that for r — > oo, V(T, r) — > V T. We thus obtain V(T, r) from the relation 

r(r, r) = -r h ( <W>> ), (14) 

where < L 2 > denotes the value of < L(0)L + (r) > for r — > oo. 

The free energy flT3|) was recently studied on 16 3 x 4 lattices for 3 and 2+1 flavor QCD 
using improved gauge and staggered fermion actions |L6], 17| . The quark masses used in 



these studies were m/T = 0.4 for 3 flavor and m u ,d/T = 0.4 and m s /T = 1 for 2+1 flavor 
case. In our analysis we use the 3 flavor potential, for which the analysis is the most 
complete. However, we have verified that differences between the potentials as functions 
of tT calculated in 2+1 and 3 flavor cases are in fact small. The resulting potential in 3 
flavor QCD (FJ) for T > T c is shown in Fig. ^| for some representative temperatures. It 



is seen that beyond a certain separation distance r = ro(T), the free energy of the heavy 
quark system becomes a constant, indicating that within the accuracy of the calculation, 
the QQ interaction potential vanishes. The resulting tq(T) as function of the temperature 
is given in Fig. ||[ Since for r > r (T), there is no more interaction between the static 
colour sources, r (T) provides a natural limit to all bound state radii. 

The potential V(T, r) obtained from Eq. (|14D is the average over color singlet and 
color octet contributions; it can be written in the form 

V(T,r) = -T In jiexph^rVT] + exp[-V 8 (T, r)/T]| , (15) 

where V\_(T,r) and Vs(T,r) specify the singlet and octet contributions, respectively. In 
perturbation theory, the leading terms for both are at high temperature and small r 
(r << T _1 ) of Coulombic form, 

K(T,r) = -±222, K 8 (T, r) = + i^, (16) 
3 r or 

with a(T) for the temperature-dependent running coupling. In the region just above the 
deconfinement point T = T c , there will certainly be significant non-perturbative effects 
of unknown form. We therefore first consider the high temperature regime, which we 
somewhat arbitrarily define as T > 1.45 T c . In this region, we attempt to parameterize 
the existing non-perturbative effects through a conventional screening form, replacing Eq. 
©by 

- -V.iT.r) = 6 V 8 (T,r) = ^ exp{-//(T)r}, (17) 
4 r 
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where /x(T) denotes the effective screening mass in the deconfined medium.^ In Fig. [| we 
show the fits obtained with the form (p~5|/|TTD , assuming a(T) and /x(T) to be unknown 
functions of T. The functional form of the potential is seen to be reproduced very well, 
with the values of a(T) and fi(T) as given in Figs. [| and |6|. We conclude that in the 
high temperature regime, a perturbative description modified by color screening gives an 
excellent account for the lattice results. The screening mass becomes a constant in units 
of the temperature, 

/i(T) = (1.15 ± 0.02) T. (18) 
A similar behavior of the screening mass was found in pure SU(2) and SU(3) gauge theory 



24| - |26] 



We shall now assume that the form fll~8 ) of the screening mass continues to remain 



valid as we lower the temperature to T c . Such a constant screening mass down to T c is 
again expected from studies of pure gauge theory. There the screening mass determined 
from the color averaged potential decreases for T < 1.5 T c as T — > T c |27||; however, the 



screening mass determined from the color singlet potential appears to be temperature 
independent for T < 10 T c H - ||. 



On the other hand, quenched QCD (pure SU(3) gauge theory) studies indicate that 
when T is lowered to T c , the perturbative ratio V\/V% = —8 will increase in favor of 
the singlet potential P%| . We therefore try to describe the behavior just above T c by a 
potential of the form fllS]), in which the color octet potential is given by 

V 8 (T,r) = ^^exp{-/,r} (19) 

instead of Eq. (O); the factor c(T) < 1 accounts for the expected reduction of octet 
interactions as T — ► T c . In the interval T c < T < 1.45 T c we thus fit the lattice results 
for V(T,r) in terms of the two parameters c(T) and a(T), with fJt(T) given by Eq. (|18|). 
In Fig. [7] it is seen that this in fact leads to an excellent parameterization of the lattice 
results ||16|| . The resulting behavior of c(T) is shown in Fig. || that of oc(T) is included in 
Fig. ||. From these considerations, we obtain the form ( jT7|) for the color singlet potential, 
in which /i(T) is given by Eq. (|i~8D and a(T) by Fig. |5j. 

4. Quarkonium Dissociation in Hot Media 

In the absence of any medium, the masses and radii of the different charmonium and 
bottonium states are quite well described by non-relativistic potential theory 0-0, based 
on the Schrodinger equation 



2m a + — V 2 + V 1 (r) 



= M?$l (20) 



where a = c,b specifies charm or bottom quarks, i denotes the quarkonium state in 
question, and r is the separation of the two heavy quarks. Above the deconfinement 



1 The form ( |l7| ) is actually only valid for r ^> 1/T. For the distances r < 1/T most relevant for 
quarkonium studies, the screening is determined by the momentum-dependent self-energy Hqq(uj = 0,p) 
p2[ . In coordinate space this leads to an r-dependent effective screening mass p3| ]. In our case this 
r-dependence of the effective screening mass is not important because of the insufficient accuracy of the 
lattice results, which in 3 flavor QCD can be fitted well with a constant effective screening mass. 
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point, the string tension vanishes and we are left with the singlet potential Vi(T,r) as 
determined in the previous section. We therefore calculate the bound state radii of the 
different states, using Eq. ( |17D in the Schrodinger equation. The radii thus obtained 
then have to be compared to the limiting binding radius ro(T) in order to specify the 
dissociation points of the different states. 

Before doing this, we have to specify the temperature scale T c . In the physically 
interesting case of 2+1 flavor and physical quark masses it is not very well determined. 
In the chiral limit of 2 flavor and 3 flavor QCD, one finds T c = (173 ± 16) MeV (2 flavor) 
and T c = (154 ± 16) MeV (3 flavor), where the quoted errors are the sum of statistical 
and estimated systematic discretization errors fT6j. For finite u and d quark masses, 
the critical temperatures coincide in 2+1 and 2 flavor QCD. We therefore assume that 
T c = 173 Mev is the relevant critical temperature and use it in our calculations. We have 
also studied the results obtained with the three-flavor temperature T = 154 Mev; the 
difference turned out to be negligable within the present accuracy of our approach. 

In Fig. |9], we illustrate the determination of the dissociation points by comparing the 
radii of the charmonium and bottonium ground states with the limiting binding radius 
ro(T). It is seen that the J/ip is dissociated at T ~ 1.1 T c , the T at T ~ 2.5 T c . These 
results are remarkably consistent with those obtained previously [13] using a screened 
Cornell potential together with lattice estimates for the screening mass. 

We now want to compare the various thresholds for quarkonium dissociation. In Table 
3, we list the different charmonium states, together with the corresponding dissociation 
temperatures. For the J/ip, this is the value T ~ 1.1 T c determined above. The \ c and 
ip' radii exceed ro(T c ), so that these two bound states cannot exist for T > T c ; the same 
holds for the x?>(2-P) and the T(3S') states. On the other hand, the dissociation points of 
the x&(l-P) an d the T(2S) states coincide approximately with that of the J/ip; all three 
values are here found to be slightly above T c . Bearing in mind the unknown systematic 
errors of the method used here, based on the intersection of r (T) and the bound state 
mass calculated with a screened Coulombic potential, it seems possible only to conclude 
that these states are dissociated very close to the deconfinement point T c . This is also 
consistent with the conclusions reached in the study of the dissociation pattern for T — > T c 
from below. The T, however, clearly persists up to temperatures well above T c . In Table 
3, we also list the screening masses relative to the deconfinement value /i(T c ) for those 
states which survive beyond T c and are then dissociated by color screening; this provides 
an indication of the change in the effective screening radius. 

Let us next address briefly the survival of those states (ip', Xc T", x'b) which cannot 
exist for T > T c . In earlier studies ]12], [TJ|], it had been assumed that they are dissociated 
at T = T c . However, a recent analysis |13|], based on the same lattice study as used 
here | IB]| , shows that they will in fact decay into open charm/beauty states already below 



T c . Such a decay becomes possible because the open charm/beauty threshold in a hot 
confining medium decreases with temperature faster than the masses of the corresponding 
quarkonium states. The relevant decay temperatures were determined in |n| and are 
included in Table 3. 

Using the results of Table 3 together with the feed-down fractions of Tables 1 and 2, we 
obtain the suppression patterns shown in Figs. pi] and llj. The new analysis of quarkonia 
in confined media has thus led to a modified three-step pattern for J/ if) suppression; only 
the suppression of directly produced J/ip requires the onset of deconfinement. For the T, 
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we obtain the multi-step form shown in Fig. 11. Here again, two states (T", x'b) decay 
below T c ; the next two, T' and Xb, are dissociated at or just above T c , and only the 
T survives much further. 

Lattice studies of the kind used here to obtain the heavy quark potential can in 
principle also determine directly the energy density e of the medium at each temperature, 
so that we should be able to give the dissociation points as well in terms of e. At present, 
however, various uncertainties (the precision of the T c determination, the quark mass 
dependence of T c and of e, dependence of e on the number of flavours, finite lattice 
effects) lead to uncertainties of about a factor two in e. A precise determination of this 
quantity is thus evidently one of the main reasons for insisting on increased computer 
performance in finite temperature lattice QCD studies; to illustrate: a 10 % error in T c 
leads to a 50 % error in e(T c ). 



qq 


T/T c 




J/tf 


1.10 


0.91 


Xc(lP) 


0.74 


* 


$(2S) 


0.1-0.2 


* 


T{IS) 


2.31 


0.43 




1.13 


0.88 


T(2S) 


1.10 


0.91 




0.83 


* 


T(35) 


0.75 


* 



Table 3: The dissociation parameters of different quarkonium states, as obtained by 
color screening for T > T c (present work) and through decay into open charm/beauty for 



Conclusions 

In the present work, we have studied quarkonium dissociation by color screening in 
a deconfining medium. Using new lattice results on the color-averaged potential in full 
QCD, we have determined the temperature dependence of the color singlet potential. 
Solving the Schroedinger equation for heavy QQ bound state with the extracted color 
singlet potential, we have specified the temperatures at which different bound states 
dissolve. Combining the results of this analysis with those recently found for quarkonium 
dissociation in a confining medium leads to the suppression patterns summarized in Figs. 



10 and 11 



For a more accurate determination of the quarkonium suppression patterns, it would 
be desirable to carry out direct lattice studies of the color singlet potential and of its quark 
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mass dependence, which may become important near the critical temperature. Further- 
more, to make contact with nuclear collision experiments, a more precise determination 
of the energy density via lattice simulations is clearly needed, as is a clarification of the 
role of a finite baryochemical potential. For the latter problem, lattice studies are so far 
very difficult; nevertheless, a recent new approach p9f could make such studies feasible. 



Finally we note that in hot media the interaction of a quarkonium state with partonic 
constituents, in particular gluon scattering, can obviously also lead to its dissociation 



30| , |31fl . A study of sequential suppression in such a framework would certainly be of 



considerable interest. 
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Appendix 



The NRQCD factorization formula for the inclusive T(nS) cross section states pT| : 



a[T(nS)} inc = a[bb 1 ( 3 S 1 )}(0 1 ( 3 S 1 ))^ +J2^bi( 3 Pj)V 1 



m 2 b 



+a[bb s ( 3 S 1 )}(0 8 ( 3 S 1 ))^ S) + a[bbsCS )](O 8 CS ))^ S) 
+ (£(2J + lyibbsCPj)}) ( ° 8(3 ^- 5) , (21) 



mi 



where <t[&6i i8 ( 2s+1 Lj)] are the sort distance cross sections, (Oi ;8 ( 2s+1 Lj))- 1 ^, nS ' ) are the 
inclusive NRQCD matrix elements, which are non-perturbative |21|, and mj, is the mass of 
the b-quark. The indices 1 and 8 refers to color singlet and color octet states, respectively. 
In the present analysis the short distance cross sections are not calculated in perturbation 
theory, but extracted from the experimental data on inclusive cross sections &(Y(nS)) inc 
flOfl , using the information on matrix elements presented in |21], [32]. The color octet 



matrix elements for J = states, i.e. (0 8 ( 1 S'o)) m [. r, " S ^ (O 8 ( 3 P )) m | nS, ' > , are compatible with 
zero within errors |^TJ. In fact it was shown that a good description of the experimental 
data on bottomonium production can be obtained simply setting these matrix elements 
to zero pUfl. Furthermore (Oi( 3 Po)) m c™ 5 ' ) is very small and we set it to zero in what 



follows. We also assume that a[bbi( 3 P 2 )]/a[bbi( Pi)] = 8.3, as predicted by perturbative 



calculations |53| p4| . With these assumptions, the inclusive cross section for S-states can 
be written as 

a[T(nS)] inc = cr[66 1 ( 3 S 1 )](0 1 ( 3 5 1 )) 1 T i n5) + cr[66 8 ( 3 5 1 )](0 8 ( 3 5 1 )) 1 T n i n5) + 

^(S.S^ 3 ^))^ + (OifPx))^), (22) 
m b 

with up = ct[66 1 ( 3 P 1 )]. The inclusive color singlet matrix elements are defined as 
(Oi( 3 Si)>; c nS) = E (0 1 ( 3 5 1 )) T(n5) P(T(mS) - T(nS)) 

m>n 
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(Oi( 3 P,)) i 



T(raS) 
inc 



E (Oi( 3 Pj))»' ( ^0MmP) - T(nS)), 



(23) 



where B(H — > T(nS)) are the inclusive branching fractions [2~I|, and (Oi( 3 S'i)) T ^ n5 - ) and 
(0i( 3 Pj))»j( mP ) are the direct color singlet matrix elements. The latter can be related 



to the bb wave function (or its derivative) at the origin [20]; we have taken this from 
|| , calculated for Buchmiiller-Tye potential. The inclusive color octet matrix elements 



(0 8 ( 3 Si)) in c were taken from f3~2| . Given the inclusive matrix and the Tevatron data 



on IS, 2S and 3S inclusive production cross sections, one can extract the short distance 
cross sections cr[66i( 3 Si)], a[bb$( 3 Si)] an d op. From the short distance cross sections we 
can estimate the direct cross section for different quarkonium states, 



(7 d [T(nS)) = 
Vd[Xbj(nP)] 



a[bb 1 ( 3 S 1 ))(0 1 ( 3 S 1 )) r ^ + a[bb 8 ( 3 S 1 )](0 8 ( 3 S 1 )) 



T{nS) 



+ 



a[bb^Pj)]{0, 



( 3 Fj)) Xw(nP) /m^. 



(24) 



In the direct cross section for the Xbj states we have neglected the color octet contribution 
which is proportional to (0 8 ( 3 Si)) Xi>j( - nP ' ) , in accordance with the results of Q32]] , where it 
was show that Xbj states are dominantly produced by a color singlet mechanism. Also the 

^ Sl ^ Xb j(nP) 



analysis of Ref. ]^TJ shows that the matrix elements (Os( 3 Si)) Xbj(nIJ) ~ 0. To complete 
our analysis we have to estimate the direct color octet matrix elements entering in Eq. 
([24]). With the assumption that (08( 3 Si)) XbJ ^ nP ^ = we can write for the direct color 
octet matrix elements 



T(15) 
inc 



(OsCS,))^ + B(T(25) - T(15))(0 8 ( 3 5 1 )) T ( 25 ) + 
B(T(3S)^T(1S))(0 8 ( 3 S 1 )) 



T(35) 



,T(35) 
' inc 



(O s ( 3 S 1 )) i 



(0 8 ( 3 ^)) T(25) 
(0 8 ( 3 5 1 )) T(35) 



B{T{3S) -> T(25))(0 8 ( 3 5 1 )) T(35) + 



(25) 



Finally, multiplying the direct cross section for different bottomonium states by the corre- 
sponding inclusive branching fractions from [2I|, we obtain the feed-down fractions given 
in Table 2. 
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Figures 
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Figure 1: The transverse momentum dependence of the inclusive production cross sections 
for different (nS) bottomonium states; the lines are exponential fits. 
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Figure 2: The color averaged potential at T = 1.03 T c (triangles), T = 1.45 T c (circles) 
and T = 3.84 T c (squares). The vertical lines indicate the points r beyond which the 
potential becomes r-independent. 
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Figure 3: The temperature dependence of the points tq beyond which the potential be- 
comes r-independent. The solid line represents the fit to the data on tq(T). The dashed 
lines represent the error band and were estimated by fitting the data on r shifted up 
(down) by one standard deviation. 
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Figure 4: The color averaged potential for temperatures T > 1.45 T c ; the lines show the 
fits described in the text. 
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Figure 5: The temperature dependence of the coupling constant; the line is a fit using the 
1-loop running coupling constant formula. 
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Figure 6: The temperature dependence of the screening mass; the line shows the average 
value. 




Figure 7: The color averaged potential for low temperatures T c < T < 1.45 T c ; the lines 
show the fits described in the text. 
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Figure 10: The J ftp suppresssion pattern. 
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Figure 11: The T suppresssion pattern. 
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